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ABSTRACT 

We determined the age of the stellar content of the Galactic halo by considering main-sequence turn-off stars. From the large number 
of halo stars provided by Sloan Digital Sky Survey, we could accurately detect the turn-off as a function of metallicity of the youngest 
dominant population, which was done by looking at the hottest (bluest) stars of a population. Using the turn-off temperature of 
a population of a given metallicity, we looked for the isochrones with that turn-off temperature and metallicity and found no age 
gradient as a function of metallicity. This would mean that this dominating population of the Galactic halo formed rapidly, probably 
during the collapse of the proto-Galactic gas. Moreover, we could find a significant number of stars with hotter temperatures than 
the turn-off, which might be blue horizontal branch (BHB) stars, blue stragglers, or main sequence stars that are younger than the 
dominant population and were probably formed in external galaxies and accreted later on to our Milky Way. 

Motivated by the current debate about the efficiency of gravitational settling (atomic diffusion) in the interior of old solar-type stars, 
we used isochrones with and without settling to determine the ages. When ignoring diffusion in the isodirones we obtained ages of 
14-16 Gyr. This result is a strong argument against inhibited diffusion in old halo field stars, since it results in a serious conflict with 
the age of the Universe of 13.7 Gyr. The age obtained including diffusion in the isochrones was 10-12 Gyr, which agrees with the 
absolute age of the old globular clusters in the inner halo. 
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1. Introduction 

The debate of how the Milky Way was for med is main l y base d 
on two classic works. The first one, by lEggen et"an (Il962h . 
states that the Galactic halo formed during a rapid, monolithic 
collapse of the proto -Galactic gas cloud. The second one, by 
' ISearle & ZinnI (Il978h . claims that the halo was formed via ac- 
, cretion of small galaxies over several Gyr in a rather chaotic 

■ manner These contrasting theories have been the subject of dis- 
cussions during the past few decades (see e.g. the reviews of 
iMajewski 1993; Helmi 2008). The knowledge of accurate and 

■ absolute ages for the various stellar populations plays an impor- 
' tant role in finding the timescale for the formation of the Galaxy. 

When determining the age of a star with isochrones, knowing 
' its accurate chemical composition, effective temperature, sur- 

■ face gravity, and distances is imperative for placing the star in 
' a restricted region of the Hertzsprung-Russell diagram (HRD). 

Stellar evolutionary models must also be accurate in order to 
date the star in a restricted time frame. 

For halo stars, distances are generally poorly known, and 
determination of the stellar atmosphere parameters has associ- 
ated errors as well. Uncertainties in the treatment of superadia- 
batic convection and gravitational settling of heavy elements can 
also influence the r elation between t he age of a s tar and its po- 
sition in the HRD llChabovej[T995t IWeissll2002h . Determining 
the age of individual stars is therefore also problematic. These 
ages are rarely known to an accuracy bett er than twenty per- 
cent (iGustafsson & Mizuno-Wiedneill200ll) . This is especially 



true for old, distant stars in the Galactic halo, constituting a big 
obstacle in the exploration of early Galactic history. 

Until recently, globular clusters ( GCs) were the main 
tracers of the Galact ic halo, (see e.g . Searle & ZinnI Il978l 
IChabover et"an Il996t ISaraiedini etall 1 19971) This picture 
has changed in the past few years especially thanks all 
the stars collected by current stellar surveys. Examples are 
the Geneva-Copenhagen survey of the solar neighborhood 
( No rdstrom et al.i ,2004i) . the HK objective-prism (Beers et akl 
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1993), the Hamburg/ESO surveys (Christlieb etaDl200lh . flie 



Sloan Digital Sky Survey (SDSS, lYork et al.ll2000 l), and the 
Sloan Ex tension for Galac tic Understanding and Exploration 
(SEGUE. I^nnv et alJl2009h . 

The results of the analyses of these surveys, combined with 
the higher confidence in obtaining accurate GCs parameters, 
has allowed a common conclusion to b e reached for the forma - 
tion of the Gala c tic hal o. For example, ISalaris & WeissI (l2002h . 
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(l2005h. and Marin-Fran ch et all .. „ 

GCs, and 'Schuster et al.' ('2006), Caro llo et~.l (l2007l 1201 61) and 
Nissen & Schuster (2010, and refere nces t herei n) using field 
halo stars, agree that both lEggen et all (Il962h and lSearle & ZinnI 
( 1978) are partially right, i.e, one part of the Galactic halo (inner 
halo) was formed during the initial collapse, whereas the other 
part (outer halo) was formed slowly and chaotically via accre- 
tion. However, the importance of the latter infall and accretion 
still needs to be quantified and the knowledge of absolute ages 
for the different stellar populations is still a missing piece in this 
puzzle. 

In this paper we aim to obtain ages for the field stellar 
populations of the Galactic halo. Considering that we do not 
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Fig. 1. Lower panel: Effective temperature of a SDSS field star 
population with metallicities in the range of [Fe/H] = -1.70 + 
0.15. Upper panel: 11 Gyr isochrone with metal fraction Z - 
2 X 10 which corresponds to the metallicity of the field pop- 
ulation. The turn-off" point of the isochrone is represented by the 
vertical red line, which agrees with the hottest stars of the field 
population. This agreement suggests an age 1 1 Gyr for the field 
population. 



know the distance and mass of distant halo stars with enough 
accuracy, we can estimate ages of stellar populations by looking 
at the temperature and metallicity of the main-sequence turn-off'. 

In a dominating coeval stellar population of a certain metal- 
licity, the turn-off' point (TO) coiTesponds to the stars with the 
highest eff'ective temperature and the isochrone with this turn- 
off" temperature provides the age of this stellar population, in- 
dependently of the distance and mass of these stars. The only 
requirement is to have a large sample of stars with independent 
measurements of metallicity and temperature, from where the 
hottest ones at a given metallicity can be found. 

This approach to estimating the age can be visualized in 
Fig. [1] where the upp er panel is an HRD w ith a theoretical 
GARSTEC isochrone (IWeiss & Schlattll |2008|) of 11 Gyr and 
metallicity of [Fe/H] = -1.70. The lower panel shows the tem- 
perature distribution of SDSS halo field stars whose metallici- 
ties span the range of +0.15 dex centered at the metal content of 
the isochrone. The turn-off temperature of the 1 1 Gyr isochrone 
from the upper panel agrees with the hottest stars of the field 
population. This suggests an age of 1 1 Gyr for this population. 
There is a significant number of stars hotter than this value. 
These stars can be BHB stars, stars of the same age but metal- 
poorer, blue stragglers, or younger stars of the same metallic- 
ity. They are present in all samples, contaminating the dominant 
population. It is also important to mention that an older popu- 
lation, which would have a colder turn-off" color than the domi- 
nating population, may well be hidden in this sample. Therefore, 
this method can only estimate the age of the youngest dominant 
population. 

We are aware that this procedure is affected by further un- 
certainties. The effective temperature, for example, is difficult to 
obtain more accurately than within 100 K in the models. It de- 



pends on the atmosphere model and also on internal convection 
processes, which are uncertain. An error of 100 K in the temper- 
ature can produce an error of up to 2 Gyr in the age of population 
II stars. By being aware of these uncertainties, we still can give 
an age estimate of the Milky Way halo and discuss implications 
ofor Galaxy formation. 

We first describe our sample and the temperature and metal- 
licity determination in Sect. |2] Second, we identify the main- 
sequence turn-off temperature in Sect. |3] which we use to deter- 
mine the age in Sect.H] Third, in Sect. |5] we calibrate our results 
using the absolute ages of GCs and finally our conclusions are 
given in Sect. |6] 



2. Sample and atmosphere parameters 

The stars for this study were taken from the SDSS/DR7 database 

iAbazajian et al. 2009), which also includes the SEGUE survey 
Yannv et al.l2 009). We used a sample of 100,000 stars with low- 
resolution spectra, from which the metallicity and temperature 
can be estimated. This was done by perfomiing a fitting between 
the observed spectra and a library of synthetic spectra. Because 
our sample of stars is large, a normal likelihood analysis can be 
very time-co nsuming. For that reason we used the MAx method 
(iJofre et al.ll201ft) to optimize the fitting. 

The library of synthetic spe ctra was built using ATLAS 12 
(ISbordonell20"0l ICastelUl I200I Sbordone, 2010 piiv. comm.) 
for generating the atmosphere models. They assume local ther- 
modynamical equilibrium and a plane-parallel, line-blanketed 
model structure in one dimension. The grid of models cover a 
range of -3.0 < [Fe/H] < -0.5 in steps of 0.1 dex, 5000 < 
Teff < 8000 K in steps of 50 K and 3.5 < logg < 5.0 in 
steps of 0.5 dex. The synt hetic spectra were then created with 
the SPECTRUM pac kage (iGrav & Corballvl[T994l) as described 
in iJofre et aT] (l20ltf ). with the difference that the input atmo- 
sphere models are now the ATLAS 12 ones. The reason for this 
change in our grid of models was to avoid interpolations between 
the synthetic spectra, which allows more reliable final parame- 
ter distribution. Also, this new grid has larger stepsize in logg, 
given that we were unable to detemiine log g with better accu- 
racy than 0.5 dex in Jofre et al. (2010) . Because our method of 
determining the turn-off only depends on the temperature and 
the metallicity, we decided not to attempt to constrain the sur- 
face grav ity parameter more than 0.5 dex. 

As in lJofreet all (l20ld), we considered a wavelength range 
of [3900, 5500] A for estimating the metallicity and tempera- 
ture. It contains the Call K absorption line at 3933 A, the Call 
H fine at 3963 A, and the Mgib triplet at 5183 A, which we 
used as metallicity indicators. This wavelength range also in- 
cludes the Balmer lines H/3, Hy, and H6 at 4861, 4340, and 4101 
A, respectively, which serve to estimate the temperature. The 
MAx method mainly weights the spectra concentrating the fit 
on these features to determine the atmosphere parameters. Since 
Ca and Mg are a-elements, a constant value of [a/Fe] - +0.4 
was assumed in the synthetic spectra to obtain [Fe/H] from those 
metallic lines. We are aware that this assumption might not be 
completely correct, as recently discussed by Nissen & Sc huster] 
(2010), who find an a-poor (i.e. [or/Fe] - +0.2) and an ar-rich 
(i.e. [a/Fe] - +0.4) population in the halo. This would mean 
that some of the stars in our sample could have a 0.2 dex off- 
set in metallicity, which is of similar to the accuracies obtained 
for metallicity with our method. Additionally, we do not ex- 
pect a high contamination o f the a- poor population found 
by iNissen & Schuster! (1201 Ol) . because it lies in the metal-rich 
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regime of our sample. In this metallicity region we also have 
contamination of disk stars. 

The mean eiTor for the metallicities and temperatures ob- 
tained is of 0.25 dex and 140 K, respectively. These uncertainties 
are low enough to allow the determination of ages, as discussed 
in Sect. |4] We also considered the atmosphere param eters pro- 
vided by the SEGUE Stellar Parameter Pipeline (SSPP. lLee etal] 
l2008a) in order to check the consistency of our work. 

We used color-color diagrams to select the F-G dwarfs and 
subgiants with 0.1 < (g - r)o < 0.48. This color constraint is 
based on the target selection for the FG sample from SEGUE, 
which represen ts a metalhcity unb iased random subsampling of 
subdwarfs (see lYannv et al.ll2009i for details on the SEGUE tar- 
get selections from color contstraints). The FG sample, however, 
does not include the stars with 0.1 < (g - r)o < 0.2, because a 
large number of BHB stars have this color as well. For our pur- 
pose of determining the turn-off (i.e. bluest color) of each stellar 
population, we decided to include these blue stars, even when 
a high contamination with BHB is produced. We refer to this 
sample as the G blue sample. 



3. Main-sequence turn-off 

Stars with a given metallicity value can have different tempera- 
tures distributed as seen in lower panel of FigJT] While the num- 
ber of stars increases smoothly at low temperatures, at high tem- 
peratures there is an abrupt decrease. In a dominating stellar pop- 
ulation we can interpret this drop as the main-sequence turn-off 
point (MSTO). In this section we describe the method employed 
to detect the turn-off and use metallicity-temperature diagrams 
to study the shape of the MSTO as a function of metallicity. 



3. 1 . Turn-off detection: Sobel Kernel edge-detection 

For low-mass metal-poor stars, the I mag nitude of the tip 
red giant branc h (TRGB) is almost c onstant. Lee et al. (1993"), 
iMadore & Fteedman ( 1995,) , ISakai eTal. (1996J, and Tabur et al.i 
(I2009h have used this knowledge to estimate distances of galax- 
ies. Through histograms of the luminosity function, they deter- 
mined the position of the TRGB by finding the luminosity where 
the count discontinuity is the greatest. They adopted a standard 
image-processing edge-detection technique to measure the mag- 
nitude of the tip, which is called Sobel Kernel. The main in- 
gredient is a first derivative operator that computes the rate of 
change across an edge, where the largest change corresponds to 
the edge, e.g the tip. 

In the same way that the TRGB shows an edge in the lumi- 
nosity function of galaxies, the main-sequence turn-off shows an 
edge in the temperature distribution function of stellar popula- 
tions, and the Sobel Kernel can also be applied to this case. In the 
following the edge-detect ion method is explained briefly only , 
but for further details see iLee et all (Il993h . ISakai et al.l (Il996l) . 
[Tabur et al. (2009), and references therein. 

We consider a continuous probability distribution function to 
model the histograms. The dependencies of binning in tempera- 
ture thus are avoided and the eiTors of the measurements are also 
included. The probability distribution can be defined by 
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Fig. 2. Example of Sobel Kernel technique used to detect the 
turn-off. The TO is indicated by the red vertical lines, with the 
dashed lines coiTesponding to their errors. Panel A: Temperature 
distribution of the field stars with [Fe/H] = -1.95 ± 0.15. 
Panel B: Sobel Kernel filter response given by Eq. Panel C: 
Histogram of temperature and the probabiHty function of Eq. ([T]i 
by the blue dotted line. 



where Teff.i, and cr,- are the temperatures and their uncertainties, 
respectively, and the total number of stars in the sample. The 
edge-detector filter applied to this temperature function is 



£(TefF) = (I>(Teff - ATeff) - (D(Teff + ATeff), 



(2) 



where EiT^g) is the filter response at the temperature Teft, and 
ATeff the bin size. 

Figure |2] illustrates the technique. Panel A corresponds to 
the temperatures of the stars with [Fe/H] = -1.95 + 0.15 form 
the SDSS survey. To better visualize the technique we plotted 
the temperature distribution with a histogram in panel C and 
the probability distribution with the blue dotted line. The Sobel 
Kernel response given by Eq. ^ is in panel B, where the maxi- 
mum response is interpreted as the MSTO. 

3.1 .1 . Uncertainties in tine turn-off detection 



(1) 



iMendeze t aP ('20021) and 'Tabur et all (l2009l) use the bootstrap 
method '^abu & Feigelson 1996) to estimate uncertainties in the 
magnitude of the TRGB. Since our task is similar to the TRGB 
detection, we decided to use the bootstrap method to study the 
errors as described by these authors. We consider the uncertainty 
in our detection as 3cr, which corresponds to 99% of the prob- 
ability of the star being a turn-off star. The uncertainties for our 
sample with [Fe/H] = -1.95 are plotted in Fig.|2]with dashed 
lines. 

Effect of errors in metallicity measurements 

The effec t of the me tallicity measurement accuracies of about 
0.25 dex dJofre et aLllToiO) in the final TO-detection were an- 
alyzed with Monte Carlo simulations. For a star with metal- 
licity value [Fe/H]o and error crjpe/Hjo, we gave a random 
value distributed as a Gaussian in the range of [[Fe/H]o - 
cr[pe/H]o, [Fe/H]o H- cr [Fe/H]o]- This was done for each star, and 
then we looked for the turn-off in this new resample. As for the 
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bootstrap eiTor, we repeated this process 500 times and calcu- 
lated the standard deviation of the TO-determination. 

We found a significant effect on the TO-detection errors due 
to the metallicity bin size used for the temperature distribution. 
This happens if the metallicity bin size is smaller than the av- 
eraged 2cr[pe/H]o, because a star can move from one bin to the 
other, which makes the temperature distribution change in those 
bins for that particular resampling. An effect in the TO-detection 
is produced, which is similar to the bootstrapping method ex- 
plained above. 

In the left hand panel of Fig. [3] we have plotted the er- 
ror due to temperature bootstrapping ("Boot") and the eiTors 
due to Monte Carlo ("MC") simulations of metallicity measure- 
ments. The binning in the metallicity for the temperature distri- 
bution was 0.2 dex, which is less than the 2cr[pe/H]|, accuracies of 
the metallicity measurements (0.25 dex). We can see how both 
curves behave similarly, where the eiTors in the TO-detection be- 
come larger on the metal-rich side. This happens because there 
are fewer stars (discussed below) at high metallicities, meaning 
that a resampling of them will affect the shape of the temperature 
distribution more. 

If the binning is much larger than the averaged 2cr[pe/H]|, er- 
rors of the metallicity measurements, there is less probability 
that a star is moved from one bin to the other one, and the tem- 
perature distribution at that resampling remains unchanged. This 
is shown in the middle panel of Fig. [3] where different curves 
represent the standard deviation of the TO-determination using 
different sizes in the metallicity binning. We can see how the er- 
rors become smaller with larger binning size, even at high metal- 
licities, where there are fewer stars. When the binning gets too 
large, the TO-detection becomes less accurate at the metal-poor 
border. For a last comparison we plotted the error due to temper- 
ature bootstrapping and the error due to Monte Carlo simulations 
of metallicity measurements in the right hand panel of Fig. |3] In 
this case the binning in the metallicity for the temperature distri- 
bution was of 0.8 dex, which is much higher than the 2cr[pe^H]o 
errors of the metallicity measurements. The errors due to metal- 
licity measurements become negligible when compared with the 
bootstrapping ones. We must take into account that the cut-off 
at the turn-off becomes inaccurate if the metallicity bin is too 
large. A binning size of ~ 0.5 dex is the limit where we obtain a 
precise turn-off determination and where the errors due to metal- 
licity measurements are negligible compared to those obtained 
with bootstrapping. 




bin = 0.8 dex 

bin = 0.5 dex 
bin = 0.2 dex 



Boot 

MC, [Fe/H] 
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Fig. 3. Left panel: errors in the TO-determination due to boot- 
strapping method in temperature (Boot, dashed line) and due 
to Monte Carlo simulations of metallicity measurement (MC, 
[Fe/H], solid line) considering a metallicity bin size of 0.2 dex. 
Middle panel: eiTors in the TO-determination with Monte Carlo 
simulations using different metallicity bin sizes. Right panel: as 
left panel, but considering the bin size of 0.8 dex for Monte Carlo 
simulations. 




[Fe/H] 

Fig. 4. Upper panel: TO-detection for samples with different 
number of stars. Lower panel: the bootstrap error in the TO- 
detection for each sample. The color lines correspond to differ- 
ent sizes used for the detection, with the number of stars (x 1000) 
indicated in the legend. 



3.1 .2. Uncertainties due to number of stars 

For the stars of the G blue sample, the larger their number, the 
smaller the errors become in the TO-detection. This is shown in 
Fig.m where the curves represent different numbers of randomly 
selected stars from the G blue sample. The amount of stars varies 
from 1,000 to 100,000. We plotted in the upper panel the value 
of the MSTO as a function of metallicity and the eiTors obtained 
with bootstrapping in the lower panel . The TO-value does not 
depend on the amount of stars used in the sample, but there is a 
dependency in the accuracies obtained. For the smallest sample 
the errors in the detection have the highest value. 

3.2. Metallicity - Temperature diagram 

The relation of the metallicity to the temperature of the stars can 
be confined in a metallicity - temperature ([Fe/H] - Teff ) diagram. 
Figure |5] shows an example of the diagram in the bottom panel. 



The top panel is the color-color diagram of the stars, which sat- 
isfies the color constraint described in Sect. |2] The majority of 
the stars are located in the main stellar locus in the color-color 
diagram, but we also see a second locus at blue g - r colors, 
which corresponds to the BHB stars . They are also present in 
the [Fe/H] - Teff diagram at high temperatures, but are not nu- 
merous enough to modify the Sobel Kernel edge detection. 

In the [Fe/H] - Teff diagram in the lower panel of Fig. |5] 
we represent the cut-off of the hottest stars. It is interpreted 
as the main sequence turn-off, and we determined it using the 
Sobel Kernel edge-detector technique explained above. At low 
temperatures the number of stars gradually decreases. This is 
partly due to selection effects, because our color constraint fa- 
vors dwarfs, which are located in the blue part of the main stel- 
lar locus. Another reason is evolution effects. Red giants evolve 
fast, meaning that we do not observe as many of them in the sky 
as main sequence stars. In addition, cold stars can also be early 
main-sequence stars, which are very faint and therefore more 
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Fig. 5. Top panel: color-color diagram for the G blue sample. 
Bottom panel: [Fe/H] - Teff diagram for the G blue sample. The 
turn-off temperature as a function of metallicity is plotted with 
blue color. 



difficult to observe than main-sequence stars close to the turn- 
off. 

The largest stellar population of our sample is the one at 
[Fe/H] = -1.65, coiTesp onding to th e peak of the metallicity 
distribution of halo stars (iBeers & C hristlieb 2005). At higher 
metallicities the turn-off halo stars are contaminated by disk 
stars. Towards lower metallicities, the number of stars decreases 
considerably, but they still populate the turn-off enough to em- 
ploy the Sobel Kernel edge-detector. In the whole metallicity 
domain, there is a contamination of hot stars that can be BHB 
stars, white dwarfs, blue stragglers, or younger m ain sequence 
stars (iPreston et al.lll994tKeston & Snedenl2000b . They are not 
numerous enough to produce an effect in the turn-off detection, 
but further analyses of them are important for understanding the 
halo structure and formation as a whole. 

The final values of the turn-off temperature as a function of 
metallicity are listed in Table [1] The first column represents the 



metallicities, in steps of 0.25 dex. For each step, we selected a 
metallicity window of bin size of 0.5 dex to determine the edge 
in the temperature distribution applying the Sobel Kernel edge- 
detector. At the second column of Table [1] we listed the temper- 
ature at the turn-off obtained for the MA^ parameters. The third 
column coiTesponds to the bootstrap errors in the TO detection, 
and the last two columns are the results of the age for the stel- 
lar populations at each metallicity. Isochrones with and without 
atomic diffusion were considered in the age determination (see 
Sect. m. For the errors in the age calculation, the uncertainties 
come from the bootstrapping of the turn-off detection and the 
consideration of 0.25 dex eiTor in the metallicity measurement. 
The size of the metallicity bin used is of 0.5 dex, see Fig. [3] 



3.2.1 . Comparison with SSPP 

As a first application of the [Fe/H] - Ten diagram, we com- 
pared the distributions of the parameters recovered by using 
MAy (Go fre et al. 2 010) and those obtained b y the SSPP(SSPP, 
iLee et all 2008a b: lAUende Prieto et al.ll2008l) . The comparison 
between th e M.Ax and SSPP parameters was discussed in 
iJofreet all (12010). where we showed the negligible general off- 
set in effective temperature of 61 K for stars with tempera- 
tures between 5000 and 8000 K. We also showed a consider- 
able offset of more than ~ 0.3 dex in metallicity for stars with 
-0.5 < [Fe/H] < -3.0. 

The [Fe/H] - Teff diagrams for both cases are illustrated in 
Fig. |6l where in the left hand panel we find the MA;^' parame- 
ters and at the middle panel the SSPP ones. We looked for the 
turn-off using the Sobel Kernel edge detector and compared the 
MSTO temperature as a function of metallicity. The 3cr boot- 
strap errors in the edge detection for both samples are compared 
in the right hand panel of the figure. 

The differences due to the offset in metallicity can be seen 
in the right and middle panels of Fig. |6]in the stellar distribu- 
tion along the [Fe/H] - Teff diagrams. In the metal-rich regime 
we see quite a large amount of stars for the SSPP parameters, 
while we do not for the MAx parameters. In a similar way, the 
metal-poor regime is more populated for the MA;^' parameters 
than for the SSPP ones. It is worth noticing that the offset of 
~ 0.3 dex in metallicity does not mean that every star is shifted 
by this value in metallicity. We can see in Fig. |6] how the peak 
of the distribution is slightly shifted, but the shape of the distri- 
bution is different for M.Ax and SSPP. This is seen especially at 
low metallicities, where the [Fe/H] values between MA;^ and 
SSPP pa rameters devi a te the most. We recall the recent com- 
ment by I Carollo et al.l (l2OI0l) about the low [Fe/H] estimates 
of SSPP, which are 0.2 - 0.3 dex higher than those values re- 
covered by high-resolution analyses. This would mean that the 
MA;^' parameters at the metal -poor regime yiel d consistent re- 
sults, which was also shown in IJofreet all (120 1 Oh in the test with 
high-resolution UVES spectra. 

In general, both distributions do not yield significantly dif- 
ferent values of the turn-off temperature for [Fe/H] < -1.2. 
The lack of M.Ax metal-rich stars does not allow any detec- 
tion of an edge for [Fe/H] > -1.2. In a similar way, the lack 
of SSPP metal-poor stars affects the determination of the TO 
at [Fe/H] < -2.7. The bootstrap errors become larger at low 
metallicities for the SSPP parameters and at high metallicities 
for the y\.Ax parameters. This is because the different numbers 
of stars seen in the metal-poor and metal-rich regimes for SSPP 
and MA;^ parameters. More stars yield a more accurate edge- 
detection (see Sect. l3.1.T] ). 
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172 
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15.5 ±3.22 



Table 1. Final value of the turn-off temperature given by the field halo stars of SDSS at different metallicities. 
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Fig. 6. Left panel: Temperatures and metallicities of the G blue stel lar sampl e (see Sect. |2]i for the atmosphere parameters estimated 
with MA;^'. Middle panel: As left panel, but with the adopted SSPP dLee e t al. 2008a) atmosphere parameters. Blue line corresponds 
to MSTO temperature for the MAx parameters and red line to SSPP ones. Right panel: Bootstrap errors in the TO detection as a 
function of metallicity. 



The errors of the TO-detection for the SSPP parameters de- 
crease smoothly with increasing metallicity, which agrees with 
the shape of the parameter distribution in the [Fe/H] - Tefi dia- 
gram of the middle panel in Fig.|6] where the number of stars in- 
creases gradually towards higher metallicities. The errors of the 
TO-detection for the MA;^' parameters, however, are relatively 
small for metallicities below - 1 .4, the point at which the errors 
increase abruptly towards higher [Fe/H]. This can be seen in the 
left hand panel of Fig. |6] where the stars disappear abruptly at 
high metallicities. These differences do not introduce significant 
effects when determining the ages, as is discussed in Sect.|4] 

3.2.2. Comparison with pliotometric sample 

The sequence of works called "uvby - /5 photometry of high- 
velocity and metal-poor stars" of W. Schuster, P. Nissen 
and collaborators analyzes a sample of hal o-field stars (e.g. 
ISchuster & Nissen 1988; S chuster et al.ll2006l) . This series con- 
sists of several papers about kinematics, chemi stry, and ages 
of halo -field stars. We considered the sample of ISchuster et al.l 
(12006'. hereafter called uvby sample), together with our SDSS 
G blue sample to show that the MSTO temperature as a function 
of metallicity is unbiased. 

Our stellar sample has ugriz colors and the uvby - /3 one has 
Stromgren photometry, and therefore we had to transform colors 
to join both samples. We only considered the MSTO tempera- 
tures indicated in Table [1] and transformed them to (b - y) col- 



ors using the scale of iRamfrez & Melendej ( |2005|) for dwarfs. 
Because our temperature values are very similar over the entire 
metallicity range, we use the fine grid version of the tables of 
I. Ramirez (priv. comm.), where the steps in temperature have a 
size of 50 K. 

The resulting MSTO (b - y) color as a function of metal- 
licity is plotted in Fig. |7] where we also placed the colors and 
metallicities of the uvby - p sample (Salaris, priv. comm.). The 
turn-off detected with spectroscopic temperatures in the SDSS 
survey agrees very well with the blue cut-off of the uvby - j3 
sample. 

We would like to point out the advantages of using the SDSS 
survey. The sample is huge, which means that the MSTO as a 
function of [Fe/H] is accurate. As a test, we applied the Sobel 
Kernel edge-detector (see Sect. 13. Il l to find the (b - y) color cut- 
off traced by the turn-off in the uvby - /3 sample. The values 
obtained as a function of metallicity are illustrated in Fig.|7] The 
uvby - p sample contains very few stars for some [Fe/H] val- 
ues (especially those at low metallicities) and the turn-off can- 
not be detected properly by the Sobel Kernel. For that same 
[Fe/H] value, however, the SDSS data provide a clear cut-off 
and the Sobel Kernel yields realistic TO values. 

Given the confidence obtained in the studies about field halo 
stars from the uvby - /3 sample, it is encouraging to obtain this 
excellent agreement in the blue cut off. Considering that both 
stellar samples are independent and that the analysis of them was 
performed using completely different techniques, this agreement 
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Fig. 7. (b - y) - [Fe/H] diagram of the photometric sample of 
Schuster et al. (Salaris, priv. comm.). The red line coiTesponds 
to the MSTO detennined from the SDSS sample as indicated 
in Table [T] and transfor med to Stromgren photomet ry using the 
empirical relations of Ram irez & Melendej (l2005b . There is a 
perfect agreement between the cut-off of SDSS and the uvby -/3 
samples. The blue dashed line corresponds to the MSTO deter- 
mined applying the Sobel Kernel edge-detector to the stars of the 
uvby - p sample. 



implies that we have found the turn-off temperature as a function 
of metallicity of the Milky Way halo, independently of the use of 
the SDSS spectroscopic or uvby - j6 photometric data set. These 
implications allow us to use the MSTO values in Table[T]in the 
following parts of this work because they are unbiased. 



4. Age determination 

When determining the age of a stellar population based only 
on its metallicity and turn-off temperature, the effects of grav- 
itational settling of heavy elements (also called atomic dif- 
fusion) can p roduce a considerable difference in the results 
(IChabover et a l. 1992a b; Salaris et al. 2000). Owing to settHng, 
the content of heavy elements in the surface progressively de- 
creases during the main-sequence phase as it sinks below the 
convective envelope. Hydrogen is thus displaced out of the cen- 
ter towards the surface. The less hydrogen there is to bum in the 
core, the faster the evolution of the star along the main-sequence. 
This effect is particularly important in old metal-poor stars, be- 
cause they have shallow convective envelopes and are exposed 
to diffusion for longer times. 

Although atomic diffusion is a basic physical effect, its ef- 
fectiveness in stars is still being debated. In the solar case, 
a satisfying agreement needs to be achieved betwee n a solar 
model and the seismic sun (e.g. 'Bahcall et al.l ll995h . In low- 
mass, metal-poor globular cluster stars there is some spectro- 
scopi c evidence of its presence (e.g. lKoriTet al]l2006tlLind et alJ 
120081 for the case of NGC 6397), but this evidence has been 
disputed dGratton et al.ll200TI) . Besides the change in the surface 
abundances over time (which implies that the presently observed 
metallicity has not been the initial one), gravitational settling of 
heavy elements affects the surface temperature in two ways: di- 
rectly through the composition-dependent opacity in the outer- 



most layers and indirectly through the accelerated evolution of 
the stellar core, which tends to be more compact than when dif- 
fusion is not taking place. As a result, Teff is lower in stars where 
diffusive s ettling is taking place. These two effects have been 
studied by IChabover et al. (1200 1[ ) by simultaneously including 



the age-reduction effect for GCs in the isochrones and by in- 
hibiting diffusion in the outermost stellar layers to prevent iron 
depletion from the surface. The latter mechan ism could be due to 
either "turbulent diffusion" ( Kom et alj2006l) . "rotation-induced 



mixing" dPinsonneault et al" 



L2OO2I), 

While in the case of GCs, age indicat ors can be used that 
only dep end very weakly on diffusion ( Chabover et al.' 1 992at 
.Meissner & WeissI |200 6), for field stars these effects become 
crucial when Teff (or the color) at the turn-off is being used, 
as in the present case. If gravitational settling is ignored in the 
isochrones employed (i.e Bergbusch & Vandenberg 1992), the 
absolute ages obta ined for the field stars can b e up to 18 Gyr 
dSchuster et al.l[T996i) . Similarly. Unavane et al.la996.) o b tained 
ages of 15-16 Gy r for a stellar sample of ICamev et al.l (1 1994b 
usmg iGreen et al.l (11987 ^ isochrones. These ages conflict with 
the age of the Universe (13.7 Gyr, Bennett et al. 2003). Stellar 
evolutionary models have improved over the years not only be- 
cause of considering atomic diffusion, but also better handling of 
opacities and a-enhanced chemical compositions. Th is has also 
significantly decreased the absolute ages of old stars (Chaboverl 
1995; Salaris & Weiss 1997). Moreover, color transformations 
for the isochrones can also affect the absolute ages based on 
the turn-of f color. An example c an be found in the uvby — /3 se- 
ries, where ISchuster et si] (l2006l) fou nd an age of 13 Gyr for the 
halo stars using newer isochrones of iBergbusch & V andenBergI 
(2001), which ignore diffusion but use empirical color transfor- 
mations that require in some cases ad hoc adjustments to lower 
the ages based on a good match between the isochrones and the 
GC turn-off colors. 

ISalaris & WeissI (1200?) have obtained field star ages lower 
by up to 4 Gyr when using stellar model tracks that included 
diffusion. Since t hat wo rk was hampered by having few objects, 

aris & WeissI (l200l suggest using SDSS data to investigate 
this question w ith a lar ger statistical sample. In the first part 
of this project (I Jofre et al. 2010), we analyzed the SDSS stel- 
lar spectra to obtain metallicities and temperatures of halo field 
subdwarfs. This is the necessary prerequisite for age determina- 
tions. The present work builds on this with the aim of determin- 
ing the age of SDSS turn-off stars, for both cases of including or 
ignoring atomic diffusion the stellar models. 



4.1. Age reduction due to atomic diffusion 

Following the method employed by Schuster et al.l (l2006l and 
references therein) and lUnavane et alT ( 19961) we looked for the 
turn-off temperature and metallicity of the isochrones that coin- 
cide with the cut-off traced by the main sequence turn-off in the 

S Fe/H1 - Teff diagram. We considered GARSTEC isochrones 
Weiss & SchlattL2008. Z. Magic, priv. comm. 2010) with and 
without diffusion. Hereafter we refer to isochrones with diffu- 
sion as "diff-isochrones" and isochrones without diffusion as 



"nodiff-isochrones". Settling of helium is responsible for the 
major change in effective temperature at the turn-off. Diffusion 
of other heavier elements produ ce a httl e additional cha nge t o 
what is produced by helium (Wei ss et al.ll2b 00; Chabove rl2007l) . 
GARSTEC models consider full diffusion of heli um a nd hy - 
drogen using the diffusion coefficients of Thoul et al.l (1 19941) . 
with no consideration of radiative levitation or rotational mix- 
ing. Since we are assuming that most of the stars have a chemical 
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Fig. 8. Solid line: Turn-ofF temperature range as a function of 
metallicity of the Gblue sample. Color lines: turn-off temper- 
ature and metallicity of GARSTEC isochrones. Upper panel: 
nodiff-isochrones of 11, 14, 16 Gyr reproduce the MSTO tem- 
perature of the field stars. Lower panel: diff'-isochrones of 9, 
11,12 Gyr reproduce the MSTO temperature of the field. 



composition that is a-enhanced, w e used isochrones accordi ng 
to this composition, as described in lMeissrier & Weissl(l2006l) . 
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Fig. 9. Solid line: MSTO temperature range as a function of 
metallicity of the Gblue sample. Color lines: different TO- 
isochrones of BASTI (top) and (bottom). 



The turn-off temperature as a function of metallicity for the 
G blue sample and TO-isochrone^ is illustrated in Fig. [8] The 
MSTO temperature range when considering the bootstrap errors 
(see Table [TJ as a function of metallicity is also plotted. The up- 
per panel shows the nodiff-isochrones and the lower panel the 
difF-isochrones. There is a need for different ages to reproduce 
the field main sequence turn-off temperature when considering 
or ignoring diffusion. When diffusion is inactive, the 14-16 Gyr 
isochrones are those that agree with the turn-off temperature. On 
the other hand, when diffusion is active, those isochrones that 
coincide with the MSTO temperature are the 10-12 Gyr ones. 
There is a mean difference of 4 Gyr in the absolu te age due to 
diffusion, as proposed by Salaiis & Weiss| (l200lh . The results 
obtaine d for ages from the n odiff -isochrones agree bet ter with 
those of lUnavane et al.1 (Il996h and lSchuster et all (Il996h . 

4.1 .1 . Comparison with otiier isoclirones 

As a comparison with the results obtained with GARSTEC 
models, we plotted the turn-off temperature of the BaSTI 
(ICassisi et al.ll2004: Pie trinfemi et al. 2004) and the (Yonsei- 
Yale. iDemarque et al.l l2004) isochrones as a function of metal- 
licity along our [Fe/H] - Tgff diagram in Fig. |9] We chose 
BaSTI isochrones as an example of models computed without 
atomic diffusion and the Y^ as an example for models that in- 
clude diffusion. For both sets of isochrones, we also considered 
only those with a-enhanced chemical composition. When us- 
ing BaSTI models, as illustrated in the upper panel of Fig. |9] 
isochrones with more than 14 Gyr are needed to reproduce the 
turn-off temperature in the metallicity range under study. This 
is consistent with the results obtained for the GARSTEC nodiff- 
isochrones. On the other hand, the lower panel of Fig. |9]shows 
how Y^ isochrones with ages of 9-12 Gyr are those that repro- 
duce the MSTO temperature, which agrees with the results we 
obtained using diff-isochrones. 

It is important to discuss the implications of this test for 
the debate about the effectiveness of atomic diffusion mentioned 
above. If there is no settling, or if the settling is completely in- 
hibited by turbulent mechanisms, then these stars would be older 
than the Universe. This is a strong argument in favor of diffusion 
acting in population II field stars. Only isochrones with atomic 
diffusion yield realistic ages of 10-12 Gyr for the dominant pop- 
ulation of field stars in the Galactic halo. Moreover, these results 
for the ages agree with those of GC, as is shown in Sect. |5] 

4.2. The age as a function of metallicity 

The final ages are listed in Table [1] The fourth column corre- 
sponds to the values obtained with the isochrones that consider 
diffusion and the fifth column gives the results from isochrones 
without atomic diffusion. The errors considered here are only 
those connected with the errors obtained for the TO-detection 
in Sect. |3] which are indicated in Table [T] as ct-yq. We deter- 
mined the errors for ages considering the turn-off temperature 
of TefF.TO + ctto and TgfFjo - ctto- 

The ages as a function of metallicity are displayed in the 
upper panel of Fig. [10] In the full metallicity domain, the ages 
obtained with nodiff-models are larger than those obtained with 
diff-isochrones. The difference between ages obtained with and 
without diffusion as a function of metallicity is displayed in 
the lower panel of Fig. (TO] where the mean difference is about 
4 Gyr. We are aware that this value could vary when using other 



TO-temperature of isochrones at different metallicities 
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Fig. 10. Upper panel: The age as a function of metallicity. The 
red line corresponds to the ages obtained using GARSTEC 
isochrones with diffusion and the black line to those without dif- 
fusion. Lower panel: absolute difference between both ages. 



stellar evolutionary models. The upper panel of Fig. [TOl shows 
that when considering the eiTors, no real trend is found for age 
as a function of metallicity. Although the oldest ages are found 
on the metal-poor side of the metallicity range, the decrease in 
age towards higher metallicities is negligible when considering 
the eiTor bars. The reason for the errors increasing at the bo rders 
of the metallicity domain is mostly that the uncertainties in the 
MSTO detection also increase at the borders. For the age deter- 
mination, we also considered the uncertainties of 0.25 dex in the 
metallicity measurements. To do this, we looked at the turn-off 
temperature of the isochrones with [Fe/H] - [Fe/H]pop - 0.25 
and [Fe/H] = [Fe/H]pop + 0.25, where [Fe/H]pop corresponds to 
the metallicity of the population. 
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Fig. 11. Ages for the MSTO temperature as a function of metal- 
licity. Black line indicates the results obtained using SSPP pa- 
rameters (Lee et al. 2008a) and the red line those from MA;^' pa- 
rameters (Jofre et al.. .2010). The offset between SSPP and 
MAx metalUcities produce negligible effects in the final age 
determination. 



We do not find any significant gradient in age and metallicity. 
This implies that there is a dominant halo population that is co- 
eval and has a small scatter in age of less than 2 Gyr. A similar 



result was found bv lSchuster et al.l (l2006h . who also see that ages 
of the halo field stars do not have a significant age-metallicity re- 
lation. As a consistency check, we determined the age for the G 
blue sample when considering the atmosphere parameters pro- 
vided by the SSPP. This is shown in Fig. [TT] where the age as 
a function of [Fe/H] for the MSTO temperature obtained from 
SSPP parameters and from the MAx parameters is plotted. Both 
results are obtained using diff-isochrones and agree very well. 
There is a small difference of 0.5 Gyr, with no trend in age as a 
function of [Fe/H]. This agreement shows the robustness of the 
MSTO temperature provided by the dominant halo field popula- 
tion from the SDSS sui-vey. 

The question of whether an age-metallicity relation exists in 
the Galactic halo is a longstanding problem. Much more effort 
has been dedicated to finding an age-metallicity relation for the 
GCs than for the field. The general belief today from GCs is that 
the outer halo of the Galaxy was formed slowly, rather chaoti- 
cally, with the Galaxy accreting material over several Gyr, while 
the inner halo was formed rapidly via a gravitational collapse. 
An age-metallicity relation for this scenario is expected to have 
a narrow age spread at low metallicities (representing the rapid 
collapse of the inner halo) and a broad spread at high metallici- 
ties (representing the merging of younger extragalactic systems). 
The limit between narrow and br oad spreads in ages is. how- 
ever, not set (e.g. [S alads & Weiss"2002'; ' De Angeli et al.ll200l 
iMarm-Franch et al..i,2009: Potter et al. IQldti . 

Finally, we recall that inner halo stars are a-enhanced pop- 
ulation II stars, meaning that they were formed in an environ- 
ment where supernova type II dominated the interstellar medium 
enrichment, thus the format ion timescale for these stars was 
very short (less than 1 Gyr; lHelmill2b08b . The age difference 
for the population II stars of our sample of 1-2 Gyr agrees 
wit h this shor t for mation timescale scenario. As discussed in 
Gilmore et al.l (11989 ). studies of the coiTelations between stel- 
lar kinematics and metallicities show that the stars that formed 
during the dissipational collapse are those with [Fe/H] < -1.5, 
and their age range is set by the ff-enhanced chemical composi- 
tion. Our results tell us that we are mainly obsei-ving inner halo 
stars, which were formed rapid l y on a short timescale during 
the collapse. iNissen & Schusted j2010l) find a large number of 
a-poor ([a/Fe] ~ -1-0.2) stars in the Galactic halo, which are 
interpreted as members of an accreted population. These stars 
have [Fe/H] > -1.6. It would be interesting to continue to study 
the kinematic and [a/Fe] properties of these stars in our sample. 



4.3. Connection to the disk 

In addition to halo stars, the SDSS/SEGUE survey contains a 
large number of disk stars, which usually have higher metallici- 
ties than halo stars. It is interesting to study how the halo and the 
disk components are connected in the [Fe/H] - Tefj diagram, to 
see whether it is possible to detect the turn-off, and to determine 
the age of that population. 

For this study we looked at different Galactic latitudes of the 
G blue stellar sample described in Sect. |2] whose Galactic posi- 
tions are known from the SDSS database. The [Fe/H] - Teff di- 
agrams for four groups of different Galactic latitudes are illus- 
trated in Fig. [121 where the regions are indicated in each panel. 
The pararneters u sed for this study are those of the SSPP pipeline 
(iLee et al ]l2008a') because the MAx parameters are focused on 
metal-poor halo stars and not metal-rich disk stars. The cut-off 
due to the main-sequence turn-off can be estimated using the 
Sobel Kernel for stars with [Fe/H] < -1.3, which is plotted with 
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Fig. 12. [Fe/H] - Teff diagrams for four different Galactic latitudes. Color lines correspond to the TO temperature as a function of 
metallicity, with blue for 90 > \b\ > 70°, green for 70 > \b\ > 50°, red for 50 > \b\ > 30° and cyan for 30 > \b\ > 0°. 



different colors depending on the Galactic latitude (see caption 
of the figure). 

We now discuss the shape of the [Fe/H] - Tgff diagram 
as a function of latitude. Let be the number of stars with 
[Fe/H]> -1.0 and the number of stars with [Fe/H]< -1.0. 
Near the Galactic poles (90 > \b\ > 70°) there is a lack of stars 
with high metallicities, with a ratio of N+/N- - 0.052. The num- 
ber of high-metallicity stars increases with decreasing Galactic 
latitude, whereas close to the Galactic equator (30 > \b\ > 0°) we 
find more metal-rich than metal-poor stars, with N+/N- - 1.24. 
This fits with the picture that the Galactic pole contains mainly 
halo stars - e.g, metal-poor stars - and the equator contain s 
mainly disk stars - e.g, metal-rich s tars (e.g .lde Jong et al.ll2010l). 
In the metallicity distributions of lAUend e Prieto et al. (2006), 
the intersection between halo and disk is at [Fe/H] = -1.0, 
which agrees with what we observe in the [Fe/H] - Teir diagrams 
of Fig. El 

We can also see an agreement of the MSTO temperature be- 
tween the different regions. This means that the turn-off stars 
in the halo have the same age, regardless of their location in 
the Galaxy. The disk stars at higher metallicities, however, can 
have bluer colors than the halo stars, and they do not present an 
abrupt cut-off in the temperature distribution as the halo does. 
This could be interpreted as the stars in the disk not representing 
a dominant population, and the presence of hotter stars tells us 
that there are considerably more younger stars in the disk than 
in the halo. 



5. Field stars and globular clusters 

Absolute ages of GCs are known with better confidence than the 
ages of field stars. Therefore, we used a sample of GCs observed 
with the SDSS ugriz filters to consolidate our results. The sample 
is described in Sect. 15.1! We discuss obsei-vational evidence of 
a common nature between cluster and field stars relating their 
colors in Sect. 15.21 Based on this evidence we use the cluster 
ages to constrain the ages of the turn-off halo stars in Sect. 15.4! 

5.1. Globular clusters sample 

The ugriz photometry of a sample of 17 clusters (An et al.' "200i 
hereafter An08) can be found on the SDSS web pagefl The ad- 
vantage in using this photometry is that we avoi d color transfor- 
matio ns, which can in troduce systematic errors (IWeiss & SalarisI 
1999] lAn et al.ll2009l) . We excluded NGC 2419, Pal 3, Pal 14, 
Pal 4 and NGC 7006 from the sample of An08 clusters be- 
cause they are too faint for SDSS cameras to detect those stars 
with magnitudes at the turn-off and main-sequence. In addi- 
tion, we excluded the cluster M7 1 because it is more metal-rich 
([Fe/H] = -0.81) than the metallicity domain used in this work. 
The resulting 1 1 clusters are summarized in Table |2] The metal- 
licities listed in th e third column are th ose of An08. That work 
used the scale of iKraft & Ivan j (l2003h for all the clusters ex- 
cept those with an asterisk, which use that o fE arris l (IT996l) .The 



' http://www.sdss.org/dr7/products/value_added/ 
anjohnson08_clusteiphotometry.html 
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fourth column of Table |2] contains the color of the turn-off (ex- 
plained below) and the the fifth column contains the reddening 
value from Schlegel et al. ( 1998). The last column has the value 
of the r ma gnitude of the turn-off taken from the fiducial se- 
quences of iAn et alj (l2008l) . 

5.2. Observational evidence 

Without any previous assumption concerning Galaxy structure 
and formation or stellar evolution theory, we present in the fol- 
lowing pieces of observational evidence that halo cluster and 
field stars may have a similar nature. To analyze the stars com- 
paratively, we first considered their magnitudes corrected by ex- 
tinction. For field stars, the magnitudes of the G blue sample 
were already de-reddened using the E{B - V) coefficients of 
ISchlegel et al. ( 1998). For the clusters, the values of the extinc- 
tion coefficient E{B — V) are i ndicated in Tab le |2] (taken from 
An08, which are also those of ISchlegel et al.lfl998i) . To trans- 
form these coefficients to ugr fil ters, we used th e ratios of selec- 
tive extinction R,,, Rg, and Rr of'McCall' ('2004'). 

In Fig. [13] we illustrate a comparison between the bluest 
color of the globular cluster and the field stars as a function of 
metallicity. To determine the blue edge for clusters and field, we 
adapted the Sobel Kernel technique (Sect. 13.1b to find the edge 
in the color distribution. The value of the (g - r) color is indi- 
cated in the fourth column of Table|2]as {g - r)edge (not coiTected 
for reddening). The numbers in parenthesis in Table |2]represent 
the bluest point of the cluster fiducial sequence taken from the 
Tables 9-28 of An08 because for these cases we were unable 
to find a good estimate of the blue edge using the Sobel Kernel 
technique. 

In Fig. [T3] we observe how the bluest color as a function of 
metallicity of the field stars agrees well with that of the GCs 
within the errors, with metal-poor stars having a bluer turn-off 
color than metal-rich ones. The trend of the blue edge of the 
field stars corresponds to the quadratic polynomial 

ig - r)edge,fieid = 0.467 + 0.217 ■ [Fe/H] + 0.039 ■ [Fe/H]^. (3) 
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Table 2. Globular clusters properti es taken from An et"aD(l2008l) 
Metallicities adopted are those of iKraft & IvansI (l2003h , except 
those with asterisk which have metallicity adopted by iHarrisI 
(Il996h . The blue edge of the color distribution calculated us- 
ing the Sobel Kernel technique is indicated in the forth column. 
Values in parenthesis are the adopt ed turn-off col or taken from 
the fiduc ial sequence of the clu ster (I An et al.ll2008h . The redden- 
ing from Schlegel et al. ( 1998) are indicated in the fifth column 
and the r magnitude of the turn-off of the fiducial sequence is 
listed in the last column. 
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Fig. 13. The solid line coiTesponds to the blue egde as a func- 
tion of metallicity for the G blue sample of field stars and the 
dashed line the bootstrap error. Red asterisks indicate the blue 
edge o f the dereddened {g - r) color and metallicity of the clus- 
ters of lAn et all (l2008l) . Their names are indicated at the top of 
the diagram. 



It is very similar to the trend of the GCs, which is 

(8 - OedgcGC = 0.471 + 0.217 ■ [Fe/H] + 0.035 ■ [Fe/Uf 



(4) 



The coefficients of the parabola are very similar and have an 
offset of 0.04 mag at the zero point, which is less than the boot- 
strapping eiTors obtained for the field stars. This suggests that 
the age obtained for the field stars should also agree with that of 
the clusters. 



5.3. Color distribution 

We selected the stars with metallicities of ([Fe/H]Gc - 
0.15, [Fe/H]Gc + 0.15), where [Fe/H]Gc is the metallicity of the 
cluster, and calculated the probability function given by Eq. ([T) 
of the color distribution. The same probability function was de- 
termined for the color distribution of the clusters. Comparisons 
for the 1 1 GCs are displayed in Fig. [14] The clusters and the field 
population at the cluster metallicity usually do not have the same 
number of stars, therefore we scaled the distribution to compare 
them better. In Fig. [14] the metallicity and the name of the clus- 
ters are indicated in each panel. 

To understand the color distributions better, we took the clus- 
ter NGC5466 as an example and plotted it separately in Fig. [TS] 
The first similarity seen in Figs. [14] and [15] is the shape of the 
distributions. This is especially notable at the turn-off, shown in 
Fig. [15] The color of the maximum number of stars in the GC 
agrees with the respective field stellar population, as can also be 
observed in Fig. [14] There are exceptions, like NGC4147, where 
the color distribution of the cluster is bluer than for the field. 
Even in this case the shape of both distributions agrees. 

Towards the red side of the color distribution, the number of 
stars decreases in both distributions. The decreasing slope for 
the field is steeper than for the clusters. This is mainly from 
selection effects. The G blue sample is biased towards turn-off 
stars. In addition, the field halo dwarfs obse rved by SPSS usu - 
ally have magnitudes of r ~ 19-20 mag (lYannv et al.ll2009l) . 
Since the SPSS survey cannot observe stars below r ~ 22 mag 
(lAbazaiian et al.ll2009l) . the samples from this survey have only 
a few distant halo stars in the main-sequence. Figure [14] shows 
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Fig. 15. Upper panel: NGC5466 color-magnitude diagram. 
Lower panel: color distribution of the cluster (dashed line) and 
the field of cluster metallicity (solid line), smoothed with the 
probability function given by Eq.[T] 
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Fig. 14. Color distribution of the GCs of Table |2] with dashed 
line. The distribution for the field stars with the same metallicity 
as the cluster is plotted with a solid line. 

us that those clusters with a luminous turn-off (see Table |2]l have 
a well-populated region on the red side of the color distribution. 
The opposite is observed for clusters with a less luminous TO. 
Finally, clusters with strong contamination of field stars, such 
as M2, M3, and Pal5, have very noisy CMD (see the color- 
magnitude diagrams in An08) and show a very populated red 
region. 

Another feature of the red side of the distribution is a minor 
peak at (g - r)o ~ 0.45 seen in each panel of Fig. [14] These stars 
are enclosed by the red dashed box in Fig. [T5]and can be found 
in the globular clusters quite easily. We can separate the dwarfs 
from the giants of the cluster by using the fiducial TO magnitude, 
as illustrated in Fig. [16] We calculated the color distributions for 
the dwarfs and giants separately, as shown in the lower panel of 
Fig. [16] Only the giants present peaks at red colors. This peak is 
the result of the accumulation of stars at the base of the red gi- 
ant branch (Iben 1967, hereafter BGB). The high concentration 
of stars at this color is an optical effect, because the base of the 
giant branch is the point where the red giant branch is the steep- 
est and the stars have the same temperature, regardless of their 
mass and evolutionary stage. Stars that lie at redder colors than 
the BGB are more massive and evolve faster. Therefore, few of 
these stars are found in the cluster, and this second agglomer- 
ation of stars then decreases again to the red end of the color 
distribution. Only the temperature of the BGB is independent of 
mass. The luminosity varies for stars with different masses. This 
is the reason there is no accumulation of stars at the BGB in the 
luminosity function plotted on the left side of Fig. [16] 



Some contamination with BGB is expected in the G blue 
sample, because their colors lie in the color constraint of this 
sample (see Sect. |2]i. An example is shown in Fig. [17] where 
color-color diagrams of the cluster M13 and the field stars with 
M13 metallicity are illustrated in the left and right panels, re- 
spectively. The blue line coiTesponds to the fiducial sequence. 
The base of the red giant branch is indicated in red. It is possi- 
ble to see how red giant branch stars are included in the Gblue 
sample due to our color constraint (see Sect. ]2}- 

5.3.1. Consequences 

The excellent agreement of the blue edge as a function of metal- 
licity of clusters and field stars seen in the [Fe/H] -color dia- 
gram of Fig. [T3] supports the idea that we are finding the main- 
sequence turn-off for a population in the Galactic halo. This 
statement is reinforced when looking at the color distribution of 
the field and GCs. The majority of the stars are located close to 
the turn-off for both field and cluster stellar populations, and the 
rate of decrease towards blue and red colors is similar. In addi- 
tion, the minor peak seen on the red side of the color distribution 
corresponds to a sample of red giant branch stars of the same 
population. 

A consequence of these results is the possible common ori- 
gin of GC and halo field stars, although our color distributions 
do not reveal the formation scenario. In particular, we can see 
here that field and cluster stars in our study may have been 
formed at the same time and in similar chemical environments 
(iBromm & Larsonll2004 . 

Another consequence, probably more extreme, is that a 
significant number of stars in the Galactic halo were formed 
in GCs. Although simulations of the dynamical evolution 
of GCs is still not c o mpletely clear (iD Ercole et all 120081; 
iBaumgardt et all l2008t ICarretta et a"D l2010l) . it is accepted 
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Fig. 17. Left panel: (u - g) - (g - r) color-color diagram of Ml 3. Right panel: color-color diagram of t he field s t ars wi th [Fe/H] = 
-1.6, according to M13 metallicity value. Blue line illustrates the fiducial sequence of the cluster of lAn et aL I (l2008l) . Red color 
corresponds to the stars at the base of the red giant branch. Red giant branch stars in the field are expected from our color constraint. 




Fig. 16. CMP of NGC 5466. The turn-off magnitude of the fidu- 
cial sequence (lAn et al. f2008) is represented by the blue horizon- 
tal line. Lower panel: color distributions of dwarfs (dotted line) 
and giants (solid line) of the cluster, separated with the limit of 
the TO magnitude. The distribution of the giants show a sec- 
ond peak, with its maximum represented by the red vertical line. 
Dashed lines indicate the width of the peak. Left panel: luminos- 
ity function of the cluster. 



that they lose stars. Obsei-vations of extratidal structures of 
GCs (Martinez Delgado et al. 2004; Law & Majewski 2010; 
Ijordi & Grebell2010l) in the Galactic halo confirm this result. But 
the question of how many field halo stars co me from di srupted 



GCs i s still under debate (e. g. Yong et afl 2008; Bol ev et al.l 



2009: Martell & Grebell l2010t ICarretta et alT r2010). Recently, 
Jordi & Grebel (2010) detected tidal features using the same data 
set as this study. This means that this particular sample of clus- 
ters is being evaporated into the Galactic halo. Although the 



extra-tidal stars of these clusters contribute a small fraction of 
the entire field halo population, they are being removed through 
tidal stripping. This fits well with the similar color distributions 
seen here. 

The case of NGC 4147 is worth mentioning . It shows a 
well-defined color distribution but shifted to bluer colors with 
respect to the field distributio n at the NGC 4147 metallicity. The 
results of iDotter et al.l (l2010h reveal that this cluster has an age 
of 12.75 Gyr, which puts it in the frame of th e normal old in- 
ner halo clusters. However. iMeissner & WeissI ("2006) obtained 
an age of 9 Gyr for this cluster, which would explain its blue 
shift. 

Finally we want to recall our statement in Sect. l4.2l that this 
sample of field stars is mainly formed from one coeval domi- 
nating population. Globular clusters are composed of stars with 
the same age and metallicitjQ. Therefore, the observations dis- 
cussed here are consistent with this picture. Moreover, we will 
see in the next section that these clusters are coeval with the field 
stars within 1-2 Gyr. 

5.4. Globular cluster ages 

From the obsei-vational evidence and implications discussed 
above, it is expected that the 11 GCs studied here have ages 
according to the dominating halo population. We discussed in 
Sect.|4]the behavior of the absolute ages obtained from the main- 
sequence turn-off. In particular where we saw how the diffusion 
parameter can reduce the age by 4 Gyr. The age of GCs is much 
better constrained, and we can use it to cons olidate the ages 
for the field stars. We considered the a ges of lAnet all (l2009l 
hereafter An09), S alaris & WeissI (l2002i hereafter SW02), and 
Dotter et al. (2010, hereafter DIO). These are shown in Table [3] 
including the metallicity scale adopted by the authors. The sec- 
ond block of the table indicates the values of An09, the third 
block the values of SW02, and the last one those of DIO. 

Although the three works use nodiff-isochrones to determine 
the ages, the methods employed differ between them. An09 and 



^ Although GC s are composed by multiple populations (e.g. 
iGratton et al.l2004h . as tracers of the Galaxy they can be treated as sin- 
gle stellar populations. 
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Table 3. A ges of the sample of globular clusters determined in three different works: lAn et al.l (l2009l An09). ISalaris & WeissI (l2002[ 



laDie i. A ges or tne sarnple or globular clusters aeterminea in tnree airterent works: lAn et al.l (IZUU^L j 
SW02) and Potter et al] (1201 Ot DIO). For the metallicity and method adopted by the authors see text. 



DIO fit the isochrones to the color-magnitude diagram near the 
turn-off, making their results being affected by atomic diffu- 
sion. The method used by SW02 to determine the ages, on the 
other hand, consists in a combination of the Ay(TO - HB) and 
A(y - /) methods (Chaboyeret al. 1992a). This approach is af- 
fected very weakly by the assumption of gra vitational settling 
(for d etails of this age-dating method, see Meissner & WeissI 
l2006h . 

The comparison between the ages of the GCs and the field 
stars is illustrated in Fig. [18] We considered the ages of the field 
stars obtained using diff- and nodiff-isochrones (see Sect. |4]l. 
Panel A shows the comparison with the An09 results. The ab- 
solute ages of the clusters agree with our results obtained with 
the nodiff-isochrones, as expected. We can see, however, that the 
trend in age as a function of metallicity is different for the clus- 
ters and the field stars, but we cannot conclude that this is related 
to a physical difference in the chemical enrichment of field and 
cluster stars. 

In panel B of Fig. [18] we find the age and metallicity of the 
SW02 clusters. In this case, we can see a similar trend between 
the age and the metallicity for the clusters and field stars. In the 
SW02 analysis, 55 clusters were included, and there was also no 
age gradient at [Fe/H] < -1.2. Interestingly, our final absolute 
age using diff-isochrones agrees with that of the clusters. This is 
further argument in favor of using isochrones with atomic diffu- 
sion for the age determination of field population II stars. 

Panel C compares the DIO results with our owns. The ab- 
solute values of the ages are different because the isochrones 
employed by DIO do not inhibit diffusion completely. They in- 
hibit settling only in the outer regions of the radiative core. In 
the outermost region, diffusion is completely stopped. In the re- 
gion belo w, diffusion has bee n Unearly ramped from zero to the 
full effect jDotter et al.ll2008l) . This partial diffusion produces a 
IVISTO temperature that is hotter than a fully diff-isochrone an d 
colder than a nodiff-isochrone (see e.g. Chabover et al.l 12001 
We downloaded a set of isochrones from the DSEP database 
and determined the age of the field stars as described in Sect. [4] 
The results are illustrated in panel D of Fig. [18] where we also 
plotted the DIO ages for the GCs. The agreement between clus- 
ters and field ages is now excellent, as expected. 

It is important to discuss here the importance of using the 
color of the turn-off for determining the ages of GCs. Not 
only atomic diffusion affects the effective temperature at the 
turn-off significantly, but other uncertainties, such as the mix- 
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Fig. 18. Pane l A: The age as a function of [Fe/H] for the clus- 
ter sample of lAn et al] ( 20091) represented by asterisks. The red 
dashed and blue solid lines are the field ages obtained using 
isochrones with and without diffusion, respectively (see Sect.[4]i. 
Panel B: Same as panel A but with asterisks representing the 
ages and [Fe/H] of Salaris & Weiss (2002). Panel C: Same 
as panels A and B , plotting this time the ages and [Fe/H] of 
iDotteret all (1201 Ob with asterisks. Panel D: Again DIO clusters 
but compared with field ages using DSEP isochrones. 



ing length parameter, metallicity scale, and color calibrations, 
can affect the absolute ages of globular clusters by 4 Gyr as 
well (Bergbusch & Va ndenBergI 12001,) . This is one of the main 
reasons for GC ages being b etter determined using th e bright- 
ness of the turn-off (see e.g. iJVIeissner AWeiss''2006^. The ab- 
solute ages obtained for the halo field stars by Schuste r"et al.l 
(2006) of 13 Gyr agree with the usual GCs ages, even when 
their isochrones do not empl oy atomic diffusion, is because the 
color transformations of the iBergbusch & VandenBergI (1200 lb 
isochrones are adjusted specifically to match the colors of the 
old metal-poor GCs. 

In summary, we have adopted three different scales for ages 
and metallicities of 1 1 clusters to calibrate the age of the Galactic 
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halo field population. We saw how they show different results, 
where An09 yield higher ages than DIO. This is because An09 
isochrones do not consider diffusion, while DIO isochrones par- 
tially do. We saw excellent agreement with the DIO scale and the 
field ages when the same isochrones were used. The SW02 ages 
employed a different method from An09, DIO, and our own. It 
has the advantage of very weakly depending on atomic diffusion 
and ad hoc adjustments in the color transformations, making it 
a more robust calibrator of absolute ages for the Galactic halo. 
It is encouraging to obtain good agreement in the absolute scale 
given by the diff-isochrones and the GCs, especially because we 
have avoid color transformations in our procedure. This rein- 
forces our statement that the inner halo has an age of 10-12 Gyr, 
where field and cluster stars were formed at the same time. 

6. Summary and concluding remarks 

We determined the age of the dominant Galactic halo popula- 
tion by using a sample of field s tars from the Sloan Digital Sky 
Survey rSDSS. lYork e"taDl2000h . 

The first main task was to determine the stellar atmosphere 
paramet ers from the SDS S sample. We used an efficient method 
(MAx, Jofre et aP l20I0h to obtain accurate parameters in a 
timely manner. Based on these parameters, we determined the 
turn-off temperature for the Galactic halo as a function of metal- 
licity, from which we could estimate the age of the halo stars. 
The turn-off was identified using a robust edge-detector kernel. 
We studied the final value of the turn-off when using the atmo- 
sphere p arameters recover ed by the SEGUE Stellar Parameter 
Pipeline jLee et alj|2008al) . We showed that it agreed with the 
value obtained using our own atmosphere parameters. In addi- 
tion, we could see the excellent agreement of the turn-off color 
obtained for the SDSS spectroscopic sample with that of the 
[Schuster et al. (2006) photometric sample. This implies that the 
edge produced by the temperature of the main sequence turn-off 
is a definite feature of the Milky Way halo, independent of the 
data set and method employed for finding it. 

For the age determination, we used the temperature of the 
turn-off and the metallicity of the stellar populations to find the 
isochrones of that particular TO temperature and metallicity. We 
explored the effect of atomic diffusion in the resulting ages and 
found an absolute difference of 4 Gyr for ages obtained ignor- 
ing or including atomic diffusion in the stellar models. The age 
was additio nally determined us ing (Demaraue et al. 2004) 
and BaSTI dCassisi et al. I I2OOI isochrones, as an example of 
isochrones with and without diffusion, respectively. We could 
see that the results obtained using BaSTI models agreed better 
with those obtained from the nodiff-isochrones, while results 
agreed better with those obtained from diff-isochrones. In this 
way, we could explain why the previous ages obtain ed with this 
method dUnavane et alj|199"6t ISchuster et al.lll996l) are greater 
than our own. Given the current debate about the efficiency with 
which atomic diffusion acts in stellar interiors, our results show 
a strong argument against fully inhibited diffusion in metal-poor 
halo stars. If fully inhibited diffusion was the case, these stars 
would be older than the Universe. Based on this argument, we 
could give an age for the field stars of the Galactic halo of 10- 
12 Gyr. 

The relation between the age and the metallicity of the dif- 
ferent stellar populations does not show a gradient. Our sample 
of stars could be a representation of the inner halo, which could 
be a product of a rapid star formation scenario, probably during 
the collapse phase of the proto-Galactic cloud. Additionally, a 
breakdown in the relation at [Fe/H] = -1.6 could be related to 



a stronger contamination of younger thick disk and outer halo 
stars. 

For the consolidation of the absolute age for the Galactic 
inner halo, we considered the ages of GCs, for which the 
parameters are better determined than for field stars. We found 
observational evidence that cluster and field stars share a similar 
history and are co mposed of a dominant population. For 11 
clusters of lAn et al . (2008), which were also observed with 
the SDSS telescopes, we compared the colors to those of the 
field stars. We could see that both color distributions agreed 
especially in the turn-off color as a function of metallicity. The 
agreeme nt betw een the ages for the field and those obtained by 
Salaris & WeissI (12002.) for GCs serves as a further argument 
that the majority of inner halo stars are 10-12 Gyr old. 

An old dominating population for the inner halo suggests that 
the Milky Way halo could have formed rapidly during a col- 
lapse of the proto-Galactic gas, where star formation took 
place 10-12 Gyr ago. This scenario agrees with conclusions 
from GC studies ('Sarajedini et al."1997' 'S alaris & Weiss1l200l 
De Angeli et al. 2005; Dotter et al. 2010), where the metal-poor 
clusters are coeval. Mo reove r , the a bsolute ages obtained for the 
clusters bv iSalaris & WeissI (|2002|) are also 10-12 Gyr, which 
implies that GCs and field stars are coeval as well. 

In our [Fe/H] - Tgff diagrams it is possible to observe a sig- 
nificant number of stars with bluer colors than the color cut- 
off of the main sequence turn-off. They could be BHB stars, 
white dwarfs, blue stragglers of metal-poor stars, but in any case 
bluer than the dominant pop ulation. They were also noticed by 
lUnavane et al.l (1199 6) and Schus ter et al.l (l2006l a nd references 
therein) and are ref erred as blue metal-poor stars (iPreston et al.l 
1994; Pr eston & S neden 2000). An interesting explanation for 
them is that they were formed in small external galaxies and have 
been accreted later on to the Milky Way halo, supporting the 
cuiTent hierarchical galaxy formation scenario. These galaxies 
experience a different star formation history to the Milky Way 
and therefore can be younger than the dominating population 
of field stars. Further analysis of their kinematics and chemical 
abundances is needed to prove this scenario. These blue stars and 
the existence of a dominant population of halo stars suggest that 
the two riva l ing fo rma tion scenarios of the Galactic halo from 
lEggen et alJ(ll962h and lSearle & Zinnl(ll978h acfiially complete, 
in a composed manner, the picture of how the Milky Way might 
have formed. These scenarios combined suggest that part of the 
halo has collapsed rapidly, while the other part has been popu- 
lated through collisions and mergers between the satellite galax- 
ies and our Milky Way. 
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